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fT^ ' Abstract. In this article we present the case of HD 41004 AB, a system composed of a KQV star and a 3.7- 

magnitude fainter M-dwarf companion. We have obtained 86 CORALIE spectra of this system with the goal 
04 ' of obtaining precise radial- velocity measurements. Since HD 41004 A and B are separated by only 0.5", in every 

spectrum taken for the radial- velocity measurement, we are observing the blended spectra of the two stars. 
^— ■ An analysis of the measurements has revealed a velocity variation with an amplitude of about 50ms~^ and a 

periodicity of 1.3 days. This radial-velocity signal is consistent with the expected variation induced by the presence 
of a companion to either HD 41004 A or HD 41004 B, or to some other effect due to e.g. activity related phenomena. 
• In particular, such a small velocity amplitude could be the signature of the presence of a very low mass giant 

Oh' planetary companion to HD 41004 A, whose light dominates the spectra. The radial-velocity measurements were 

then complemented with a photometric campaign and with the analysis of the bisector of the CORALIE Cross- 
5_( ' Correlation Function (CCF). While the former revealed no significant variations within the observational precision 

^ I of ~0. 003-0. 004 mag (except for an observed flare event), the bisector analysis showed that the line proflles are 

. varying in phase with the radial- velocity. This latter result, complemented with a series of simulations, has shown 

that we can explain the observations by considering that HD 41004 B has a brown-dwarf companion orbiting with 
the observed 1.3-day period. As the spectrum of the fainter HD 41004 B "moves" relative to the one of HD 41004 A 
(with an amplitude of a few kms~^), the relative position of the spectral lines of the two spectra changes, thus 
5— i ' changing the blended line-profiles. This variation is large enough to explain the observed radial- velocity and 

bisector variations, and is compatible with the absence of any photometric signal. If conflrmed, this detection 
represents the flrst discovery of a brown dwarf in a very short period (1.3-day) orbit around an M dwarf. Finally, 
this case should be taken as a serious warning about the importance of analyzing the bisector when looking for 
planets using radial- velocity techniques. 

Key words, techniques: radial velocities - binaries: visual - binaries: spectroscopic - stars: brown dwarfs - stars: 
exoplanets - stars: individual: HD 41004 
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1. Introduction 



RadialjVtiloclty tiicliiiiqueH hnvn m f&Y mwailad about 80 
planetiiy (!ompaiiioiiH ai'ouiid Holai' iypa dwai'feQ. The 
most precise instruments currently available for planet 
searches can measure the velocity of a star in the direc- 



tion of the line-of-sight with a precision of the order of 2- 
3ms~^ (e.g. pueloz ct al. 200 la|; [Butler et al. 2001|; |Pepe 



et al. 2002), but even higher precision is expected from 



Send offprint requests to: Nuno C. Santos, e-mail: 
Nuno . SantosSobs . unige . ch 

* Based on observations collected at the La Silla Observatory, 
ESO (Chile), with the CORALIE spectrograph at the 1.2- 
m Euler Swiss telescope and with the Stromgren Automatic 
Telescope (SAT). 

^ See e.g. obswww. unige. ch/~udry/planet/planet. html 



instruments available in the near future (e.g. HARPS - 
Pepe et al. 2000a). This will definitely allow the discov- 



ery of lower mass and longer period planets, that remained 
undetected up to now due to the low amplitude of the in- 
duced radial-velocity variation. 

The gain in precision will, however, bring to light some 
of the limitations of the radial-velocity method. It is well 
known, for example, that radial-velocity "jitter" with am- 
plitudes up to a few tens of ms~^ is expected to result from 
the presence of strong photospheric features like spots 
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Table 1. Stellar parameters for HD 41004 A 



Parameter 



Value 



Reference 



Spectral type 

Parallax [mas] 
Distance [pc] 

B-V 

Teff- [K] 

log 3 [cgs] 
Mv 

Luminosity [Lq] 
Mass [Mq] 

l0g7?HK 

Age [Gyr] 
Prot [days] 

V sini [kms""*"] 
[Fe/H] 



K1V/K2V 

23.24 ± 1.02 
43 
8.65 
0.887 
5010 
4.42 
5.48 
0.65 
~0.7 

-4.66 
1.6 
~27 

1.22 
-0.09/+0.10 



Hipparcos (|ESA 19971) / 

uvby (see Sect 2.4 ) 
Hipparcos ( lESA 19971) 
Hipparcos ( ESA 1997 ) 
Hipparcos ( ESA 1997 ) 
Hipparcos ( ESA 1997 ) 

See text 

See text 



Flower (1996) 



Henry et al. ( 1996| ) 



Donahue 



T99IT 



Noyes et al. (1984) 



CORALIE 
uvby/CORALlE 




6704 



6706 6708 

A (A) 



6710 



6712 




3964 3966 3968 3970 3972 

X (A) 



or convective inhomogeneities, associated with chromo- 



spheric activity phenomena (3aar & Donahue 1997; Saar 
et al. 1 |99^ ; [Santos et al. 200C ). The presence of spots can 



even induce a periodic radial-velocity signal similar to the 
one expected from the presence of a planet. This is the case 
for HD 166435 ( Qucloz et al. 2001b| ), a star presenting a 
radial-velocity signal with a period of about 3.8-days, but 
showing both photometric and bisector variations with the 
same periodicity. 

The case of HD 166435 illustrates very well the need 
to confirm, at least for the shortest period cases, that the 
radial- velocity signature is indeed due to the presence of a 
low mass companion, and not due to some kind of intrinsic 
phenomena. As shown by [Queloz et al. (2001h| ), the use 
of photometric data and bisector analysis was crucial to 
clarify the origin of the radial-velocity variations observed 
on this star. 

In this paper we present the case of HD 41004, a visual 
double system consisting of a K1V-M2V pair (A and B 
components). This system was found to present a radial- 
velocity signature similar to the one expected as if the 
Kl dwarf had a very low mass planetary companion in a 
1.3-day period orbit. Although the photometric data re- 
vealed no significant photometric variations, an analysis of 
the Bisector Inverse Slope ( BIS) of the CORAL IE Cross- 
Correlation Function (CCF) ( Qucloz ct al. 200l"b| ) revealed 
a periodic variation in phase with the radial- velocity sig- 
nal. In the following sections we will show that the radial- 
velocity variation is in fact not a result of the periodic 
motion of the A component, but of the Doppler motion 
of the spectrum of the B component due to the presence 
of a brown-dwarf companion. The results strongly caution 
about the need to use methods capable of detecting line 
asymmetries, like the bisector analysis, when dealing with 
high-precision planet searches with radial-velocity tech- 
niques. 



Fig. 1. Upper panel: Added CORALIE spectrum in the 
Lithium line region for HD 41004. The plot reveals no clear 
Li detection for this star. Lower panel: Added CORALIE 
spectrum in the Can H line central region for HD 41004. 
The strong emission near line center suggests that the star 
is chromospherically active 



2. The case of HD 41004 

2.1. Stellar characteristics 

In the Hipparcos catalogue HD 41004 (HIP 28393, 
CD —48 2083) is presented as the A component of a double 
system, where the 3.68 mag fainter companion is at a sep- 
aration of about 0.5". From the magnitude difference, and 
considering that both stars are at the same distance (i.e. 
that they constitute a physical double system) we can de- 
duce that HD 41004 B is probably a M2 dwarf. In the rest 
of this paper we will denote by HD 41004 AB the system 
composed by HD 41004 A and HD 41004 B. 

The small angular separation between the two com- 
ponents implies that we cannot separate (with the in- 
strumentation used) the light coming from the two stars. 
However, given the relatively "large" magnitude differ- 
ence, we do not expect the flux coming from HD 41004 B 
(about 3% of that coming from the A component) to affect 
significantly the determination of the stellar parameters of 
HD 41004 A. All the parameters described below are thus, 
in a good approximation, those of HD 41004 A (the star 
whose light dominates). 

At a distance of ^43 parsec - 7r=23.24±1.02mas - 



ESA (19971 ), HD 41004 A is a Kl dwarf shining with a 
visual magnitude V=8.65 in the southern constellation 
Pictoris (the Painter's Easel). Its colour index is B — V = 
0.887, as listed by the Hipparcos catalogue and its abso- 
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Table 2. Elements of the fitted orbit 



p 


1.3298 


± 0.0002 


d 


ai sini 


0.0009 


± 0.0000 


Gm 


T 


2452200.43 


± 0.13 


d 


et 


0.07 


± 0.04 




Vr 


42.532 


± 0.002 


kms 


a;+ 


23 


± 35 


degr 




50 


± 2 


ms~ 


Mm) 


0.1741 • 10"^° 


± 0.0224 ■ 10"^° 


M© 


a{0 ~ C) 


14 




ms 


N 


86 








t Cons istent w ith a circular orbit according to the Lucy & 
Sweene|r("l97l[) test. 



lute magnitude Mv=5.48. The basic stellar parameters are 
summarized in Table |l|. 

There are not many references in the literature on 
HD 41004 A, and no precise spectroscopic analysis is avail- 
able for this star. Using the b — y colour from Olsen 
(1994b)0 and the calibration of [Alonso et al. (1996| ) we 
obtain a temperature of 5030 K, 100 K higher than the 
4930 K obtained from the same colour but using the cali- 
bration of Olsen (1984). Using the B — V colour wc would 
obtain Teff=5075K from the calibration of Flower (1996| ). 
We have decided to use an average value in the following 
of this paper (Teff=5010K). Photometric calibrations can 
also be used to estimate the surface gravity of a dwarf. 
Using the calibration of plscn (1984 ) we have obtained 
log 5=4.54, slightly higher than (but compatible with) the 
value of 4.30 obtained from the Hipparcos parallax (e.g. 
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Fig. 2. Upper panel: Phased CORALIE radial- 
velocity measurements (subtracted by the 7-velocity) for 
HD 41004 AB and best Keplerian fit. Lower panel: O - C 
residuals of the CORALIE radial- velocity data over time. 
The lines represent the best linear fit without consider- 
ing the first measurement (solid line) and a parabolic fit 
(dashed line). The linear fit has a slope of 0.2 ms^^ day~^ 



Allendp~Pi'ic;Lu (jL al. 1999 ) - again, au average value was 
considc^Tjd. The ubLaiued Tcff and \ ogg are cu niyaLible 
with a spectral classification of KIV ( ESA 1997 ) 



firmed by the value of log i?^j^=— 4.66 obtained by Henry 



From the photometric calibration of 3chuster & Nissen 
(1989|)]we have obtained a value of [Fe/H]=-0.09. Another 
estimation of the metallicity can be obtained from the 
analysis of the Cross-Correlation Function (CCF ) of" 
CORALIE. In fact, its equivalent width (or surface) is 
very well correlated with the metallicity of the star and 
its colour B — V (Sect. |A^ ). Using this method we have 
derived a higher value of [Fe/H]=-|-0.10. 

Using a metallicity of +0.10, the temperature of 
5010 K, and its absolute magnitude, we can estimate a 
mass of ^0.7 fo r HD 41004 A using the isochrones of 



et al. (1996), and by the fact that a flare event was ob- 
served in photometry (s ee Sect.[2.4|). From t his value, and 
using the calibration of Noyes et al. (1984), we can esti- 
mate the rotational period to be ~27 days, a value per- 
fectly compatible with the low v sin i measured. Using the 
calibration of Donahue (1993D (also presented in Henry et 
al. 199q ), based on the chromospheric activity level we 



Schallcr et al. (1992] ) 



The CCF can also be used to determine the projected 
rotational velocity v smi of a star. From the calibration 
presented in Sect. A.l, and taking the average Gaussian 
width of the CCF for our star (4.36kms~^), we have 
obtained a value of v sini=1.22kms^^, similar to the 
1.54kms-i obtained from the CORAVEL CCF using the 



calibration of |Bcnz fc Mayor (1984| ) . 

An analysis of the Can H fine core of HD 41004 A re- 
veals a strong emission, suggesting that the star is chro- 
mospherically active (see Fig.Q lower panel). This is con- 



derive an age of 1.6 Gyr for this star. 

In Fig.|l| (upper panel) we show the 6708 A Li-line re- 
gion for HD41004AB. There is some hint of the presence 
of Li, but the small feature at the position of the Li line 
has an equivalent width of less than 2 mA. Considering 
the stellar parameters of our star, this implies an upper 
limit for the Li abundance of log A^(Li)<0.0 (this was de- 
rived from a LTE spectral analysis using the line abun- 
dance code M OOG ( ^ncden 1973 ), and a grid of Kurucz 
( Kurucz 1993| ) ATLAS9 atmospheres). Given the low tem- 
perature of our star, this limit cannot be used as a strong 
age constraint; the age derived from the activity level is 
nevertheless compatible with such a low value for the Li 
abundance (e.g. Jones et al. 1999). 



A similar b — y was obtained by us - see Sect. 2.4 



2.2. Radial-velocity data 

HD 41004 AB is included in the sample observed in 
the context of the Geneva extra-solar planet search 
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programme, at La Silla Observatory (Chile), using 
the CORALIE spectrograph ( |Udry et al. 20001 ). From 
December 2000 to February 2002 we have obtained a to- 
tal of 86 precise radial- velocity measurements^ of this sys- 
tem. The average individual photon-noise precision of the 
radial- velocity measurements is 5.5 ms~^. 

Since the angular separation of the two stars is much 
smaller than the diameter of the CORALIE fiber, in each 
spectrum obtained we are seeing the blended spectra of 
HD 41004 A and HD41004B. This means that the result- 
ing radial- velocity is not the one of any of the two individ- 
ual stars, but rather it corresponds to a weighted average 
of the radial velocities of the two components (although it 
will be closer to the radial velocity of the A component, 
since the M dwarf is "only" contributing to ^3% of the 
light). 

The analysis of the radial-velocity revealed a varia- 
tion with a period of ~1.3-days, and an amplitude of 
50ms-i (see Table|). A fit to the data (Fig.|- upper 
panel) suggested that the radial-velocity variation could 
be due to the presence of a low mass (0.25 Mjup) plane- 
tary companion orbiting HD 41004 A at a distance of only 
about 0.02 A. U. (i.e. ^--^5.5 stellar radius, considering that 
HD 41004 A has a radii typical for a Kl dwarf - 0.8 R©) 
in a circular trajectory. Although this possibility was very 
tempting, as we will see below the planetary explanation 
is not the best. 

The relatively high O — C residuals of the fit (14 ms^^) 
are quite difficult to explain in the light of activity phe- 
nomena, since we do not expect them to induce such a 
high radial-velocity "jitter" in a slow rotator K dwarf 
( ^aar et al. 1998| ; [Santos et al. 2000| ). Actually, part of 
the noise seems to be coming from the existence of some 
radial-velocity drift. An analysis of the (O — C) residu- 
als to the fit reveals a clear trend (considering the last 
group of points only), with a gradient of 75ms^^yr~^ - 
(Fig. ||- lower panel). We note that a comparison with old 
CORAVEL measurements, as well as with the first (iso- 
lated) CORALIE measurement, does not corroborate this 
trend, i.e. the old measurements (see Fig. |[ lower panel) 
are not aligned with a linear fit to the residuals of the 
last group of points. This result is compatible with the 
fact that the radial-velocity signal passed through a min- 
imum. Finally, some noise might also be coming from 
the influence of the companion to HD 41004 A. As shown 



by Pepe et al. (2000b), the presence of a close compan- 
ion might introduce some noise in the computed radial- 
velocity due to variations of the flux ratio of the two stars 
within the CORALIE fiber as a function of the seeing con- 
ditions. This noise depends on the separation between the 
two objects, their magnitude difference, and their radial- 
velocity difference. For this case we do not expect, how- 
ever, a strong effect, mostly because of the small angular 
separation between the two objects (the angular separa- 
tion of the two stars is much smaller than the diameter 
of the CORALIE fiber, and thus the influence should be 
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Fig. 3. Measured radial-velocity vs. BIS for HD 41004. 
The plot shows the two variables vary in phase. The best 
linear fit is shown. The slope of the line has a value of 0.67 



always similar and fairly independent e.g. of the seeing 
conditions) . 

As mentioned above, HD 41004 B (a M2 dwarf) is lo- 
cated at an angular separation of ~0.5". At the distance of 
HD 41004 AB this corresponds to a projected separation 
of ^21 A. U., considering that the two stars form a real 
system^. This value corresponds to the minimum "real 
separation" (since we are not able to know the distance 
between the two stars in the direction of the line-of-sight) 
and implies a period of ~90yr, using a mass of O.TMq and 
0.4Mq for HD 41004 A and B, respectively. Considering a 
circular orbit, the amplitude of the variation has a value 
of 2.5kms~^, and the maximum derivative of the veloc- 
ity would be 170ms~^yr~^. This value is perfectly com- 
patible with the observations, suggesting that the radial- 
velocity trend observed has its origin in the motion of 
HD 41004 A around its lower mass companion HD 41004 B. 
But we cannot exclude that (at least part of) the observed 
trend might have a different origin, like e.g. the presence 
of a low mass companion to HD 41 004 A. 

2.3. Bisector analysis 

With such a short period "orbit" (and further knowing 
that the star is chromospherically active) , it was extremely 
important to analyze the bisector of the CCF for this 
object. Such an analysis, as presented in Qucloz et al 



(20011: ) , is in fact an excellent way of discriminating be- 



tween radial- velocity variations due to changes in the spec- 



Available in electronic form at the CDS. 



^ This is corroborated by the simulation presented below, 
that implies that both stars have similar radial-velocities. 
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Fig. 4. Magnitude differences HD 41004 AB-HD 39823 in 
the instrumental system. The horizontal lines represent 
the mean values (JD 2452276 excluded), see text 



tral line shapes from real variations due to the orbital mo- 
tion of the star. In a few words, we compute the (bisec- 
tor) velocity for 10 different levels of the Cross-Correlation 
Functiorj^. The values for the upper (near continuum) and 
lower bisector points are averaged and subtracted (see 
Fig.^. The resulting quantity, the Bisector Inverse Slope 
(BIS), can be used to measure the variations of the line 
bisector. 

The result of the bisector analysis for our star is pre- 
sented in Fig.^. As can be seen from the figure, the 
Bisector Inverse Slope (BIS - Queloz et al. 2001b) - simi- 



lar to the usual bisector span (Gray 1992) - varies in phase 



with the radial-velocity. This indicates that most proba- 
bly the radial- velocity variation observed is being induced 
by some "intrinsic" phenomena, and not by the presence 
of a planet around HD 41004 A. 

We also note that the radial-velocity in this plot was 
computed using a different cross-correlation mask than 
the one used to determine the radial- velocities pres ented 
in Fig.|| (for more details see Queloz et al. 2001b ), and 



^ We can define 12 levels in total, i.e., dividing the CCF in 
11 slices, but for clear reasons we exclude the continuum level 
and the level corresponding to the tip of the CCF, 



Fig. 5. Magnitude differences HD 41004 AB-HD 39823 in 
the instrumental system on JD 2452276, The horizontal 
lines represent the mean values (JD 2452276 excluded), 
see text 



the resulting radial- velocity variation amplitude is smaller 
than the one measured with the standard mask (37 ms"^ 
instead of 50 ms^^). This large difference would not be 
expected if the radial-velocity signal was induced by the 
presence of a planet around HD 41004 A, since in that case 
all the lines would "move" by the same amount, leav- 
ing the radial-velocity variation (but not necessarily the 
"zero" point) fairly independent of the mask used. The 
different amplitude observed may be interpreted as a hint 
that some phenomena is not affecting in the same way all 
the spectral lines observed. 



2.4. Photometry 

The observations described above made us suspect that 
similarly to HD 166435, our star should present photomet- 
ric variations compatible with the presence o f a rotatin g 
spot. According to the Hipparcos catalogue ( ESA 1997 ), 
however, HD 41004 AB is stable, with a magnitude scat- 
ter of 0,016 mag, a value that is typical for a star of its 
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Table 3. Standard uvby indices for HD41004AB and the 
comparison stars 



Star 


V 


(b-y) 


mi 


Cl 


N 


HD41004AB 


8.621 


0.524 


0.404 


0.313 


49 




±0.007 


±0.004 


±0.007 


±0.008 




HD38459 


8.493 


0.503 


0.368 


0.329 


30 




±0.004 


±0.003 


±0.006 


±0.007 




HD39823 


8.381 


0.526 


0.315 


0.312 


38 




±0.004 


±0.003 


±0.007 


±0.006 




HD43548 


8.805 


0.434 


0.225 


0.309 


26 




±0.005 


±0.004 


±0.007 


±0.006 





magnitudcQ But this hmit is not very teUing, since a spot 
with a fiUing factor of a few percent can already induce a 
radial- velocity signal with an amplitude of a few ms""'^ in 



our star (3aar & Donahue 1997). Furthermore, we could 
suppose that this star had just developed a spotted re- 
gion, that was not present by the time of the Hipparcos 
measurements . 

In order to control the photometric stability]], we have 
collected precise photometric data for HD 41004 AB. 

The photometric uvby observations of HD41004AB 
were obtained at the Stromgren Automatic Telescope 
(SAT) at ESO, La SiUa, Chile on 11 nights between 
December 2001-February 2002. Details on the spectrom- 
eter and the fully-automatic mode of the telescope are 
given by plsen (1993| , ^994^. 

Three comparison stars, HD 39823 (Cl), HD 38459 
(C2), and HD 43548 (C3) were observed together with 
HD41004AB in the sequence C1-HD41004-C2-HD41004- 
C3-HD41004-C1, allowing accurate differential magni- 
tude differences to be formed. Each observation consisted 
of three individual integrations, and the total number 
of photo-electrons counted per observation was at least 
75,000 in u, and considerably more in the other three chan- 
nels. Sky measurements were obtained at a position near 
HD41004AB at least once per sequence. A circular di- 
aphragm of 17" diameter was used throughout. Nightly 
linear extinction coefficients were determined from the 
observations of the comparison stars and other constant 
stars, and when appropriate, linear or quadratic correc- 
tions for drift during the nights were applied. 

Differential magnitudes (instrumental system) were 
formed using for each candidate observation the two near- 
est comparison star observations. All comparison star ob- 
servations were used with C2 and C3 first shifted to the 



Again, the photometry is most of all sensitive to the A 
component of the system, and we are thus basically measuring 
the eventual photometric variations of HD 41004 A. For exam- 
ple, if the HD 41004 B would change its luminosity by 5%, the 
total flux of the HD 41004 AB system would change by only 
~0.15%, i.e., about 0.002 magnitudes, a value that is within 
the observational errors. 

''' These observations were also done with the goal of search- 
ing for an eventual planetary transit. Given the short period 
the probability of transit would be particularly high. 
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Fig. 6. Phase folded magnitude differences 
HD 41004 AB-HD 39823 in the instrumental system; 
JD 2452276 observations have been excluded. The hori- 
zontal lines represent the mean values, and vertical lines 
are drawn at phase 0.25, i.e. at the phase of expected 
central transit 



level of Cl. A careful check of the comparison stars were 
performed, and during the observing period they were all 
found to be constant within 0.003-0.004 mag. (vby) and 
0.005-0.006 mag. (u), which is close to the observational 
accuracy. 

In Fig.^ we show the 152 obtained 
HD 41004 AB-HD 39823 magnitude differences in 
the instrumental system^. Typical rms errors of one mag- 
nitude difference are 0.004-0.005 (ybv) and 0.006-0.008 
(u). On one night, JD2452276, a sudden brightening of 
HD41004AB was observed, presumably due to a flare. 
Details are shown in Fig. ||. 

Mean magnitude differences HD 41004 AB - HD 39823, 
excluding the JD2452276 observations, are 0.240 ± 0.004 
(y), 0.238 ± 0.004 (6), 0.311 ± 0.005 (w), and 0.385 ± 0.008 
(m). Thus, the rms errors are close to those expected from 
the observational uncertainties. We see no sign of variabil- 
ity correlated with the 1.3 day period (Fig.^). However, 

* The individual measurements are available in electronic 
form at CDS. 
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Fig. 7. Average bisector for the KO dwarf HD 13445. This 
figure shows that the Une bisector is almost vertical for 
a star of this spectral type. The two regions denoted by 
the dotted lines (Viow and Vhigh) represent the intervals 
used to compute the Bisector Inverse Slope (BIS), defined 
as BIS — Vhigh— Viow: Vhigh/iow are simply the average 
of the velocity of th e 4 "points" in each of the intervals 
(for more details see Queloz et al. 200Th| ). The error bars 
represent the value ofCT(iy7y^(iVyPwhere a{i) is the rms 
around the mean velocity obtained for a given bisector 
level i from the N — 118 CCF's available. A and B denote 
the upper and lower bisector points 



as the JD2452276 observations show that HD41004AB is 
active, long term variability at some level can not be ex- 
cluded. 

As an extensive list of uvby standard stars were also 
observed, an accurate transformation to the standard sys- 
tem could be done. The selection of standard stars and the 
transformation to the standard system will be described 
in detail elsewhere. Standard uvby indices are given in 
Table|. For HD 41004 AB they are all within 0.035 mag 
of the mean uvby indices for the MK spectral classification 
K2V. The colour {b — y) only differs by 0.006 mag from 
the mean K2V colour index (cf. plsen 1984| ). 

The results, revealing that the star is stable within the 
instrumental precision^, support the idea that the radial- 
velocity variation observed cannot be due to the pres- 
ence of spots in the surface of HD 41004 A, as was the 
case for HD 166435. But it is interesting to mention that 
indeed we would not expect this to be the case, given 
that the rotational period of HD 41004 A (see Table|^) is 
much longer than the observed 1.3-day period in radial- 



We note also that no transits were found. 



Fig. 8. Figure illustrating the effect of changing the rel- 
ative position of the CCF's of the HD 41004 A (deeper 
CCF) and HD41004B (smaller CCF): the resulting 
blended CCF (the addition of the CCF's of the two stars) 
will change its profile, with the consequent variations in 
the observed bisector shape (dashed lines). The CCF's 
and bisectors in this figure are exaggerated for illustration 
purposes 

velocity]^. Otherwise, and given the low v sini for this 
star, we would be observing a system in a very odd con- 
figuration. In other words, HD 41004 A is not a HD 166435- 
analog. This is further supported by the fact that contrar- 
ily to HD 166435, the radial-velocity data varies in phase 
(and not anti-phase) with the BIS (see Fig. ||). 

3. Simulating the observations 

The observed correlation between the BIS and the radial- 
velocity signal, together with the absence of any signifi- 
cant photometric variation lead us to analyze the possi- 
bility that the radial-velocity and bisector changes were 
being induced by a variation in the relative position of 
the cross-correlation dips of HD 41004 A and HD 41004B. 
Since the angular separation between the two objects is 
only ~0.5", it is impossible to observe them separately 
with the 2" fiber of the CORALIE spectrograph. For each 
spectrum obtained for HD 41004 AB we have thus the ad- 
dition of the spectra arriving from the two components. 
The question was thus to know if a periodic change of the 

^" We note that the rotational period was obtained from the 
chromospheric activity level, this latter being derived from an 
analysis of the Call H and K lines. These lines, located in the 
blue part of the spectrum, are probably not affected at all by 
the blend with the spectrum of the M dwarf, and will thus give 
us information only regarding HD 41004 A. 



8 



N.C. Santos et al.: A brown dwarf around HD 41004 B 



relative position of the two spectra could induce the ob- 
served radial-velocity and bisector variations. Two main 
possibilities can be explored in this sense. 

3.1. Case 1 

First, we can consider that HD 4 1004 A has a very low 
mass companion that is responsible for the observed 
radial- velocity variation, and that HD 41004 B has a "con- 
stant" velocity. If that were the case, what would be the 
effect on the shape of the bisector (i.e. on the BIS)? 

A simple simulation, adding two Gaussian functions 
(i.e. two cross-correlation dips), the first contributing to 
only ~3% of the flux[] (representing HD 41004 B) and the 
other with --97% of the flux (representing HD 41004 A), 
has shown that changing the position of the deeper func- 
tion by 50ms~^ will reproduce the observed variation in 
the measured velocity, but leave "constant" the shape of 
the bisector. This showed us that we can exclude the pos- 
sibility that the observed radial-velocity variation with a 
period of 1.3 days is produced by a planetary companion 
to HD 41004 A. 

This simulation suggests that if we hope to simulate 
the large observed bisector shape variations we need to 
change the relative position of the two cross-correlation 
dips by (much) more than 50 ms""'^. 

3.2. Case 2 

To try to verify whether large amplitude periodic variation 
in the relative position of the two CCFs can induce the ob- 
served signal, we have done a set of simulations where two 
cross-correlation functions were added. As above, these 
two functions, one corresponding to HD 41004 A and the 
other to HD 41004 B, were weighted by the relative flux of 
the two stars (A 30 times brighter than B). These CCF's 
were considered to be Gaussian for a question of simplic- 
ityQ. This is a good approximation, since the bisector of a 
K dwarf is o bserved to be quite vertical (Fig. 0) - see also 
[Gray (1992|) . However, and as we will see below, this also 
poses some problems when trying to constrain the results. 

In each simulation, the radial-velocity of the primary 
dip (corresponding to HD 41004 A) was considered to be 
fixed, and the smaller secondary dip, corresponding to 
HD 41004 B, was taken to vary periodically with a given 
velocity amplitude {K2) (see Fig.^). This variation is sup- 
posed to simulate the presence of a companion around 
HD41004B. 

A grid of simulations was done by changing both the 
velocity amplitude K2 (in the range 2-10kms~^) and the 
parameters of the "small" dip (the width (T2 , and the depth 
A2 were changed in the range 5-10 km s~^ and 0.04-0.20, 
respectively) , as well as the difference between the velocity 




0.26 0.5 0.75 10 0.25 0.5 0.75 
Phase Phase 

Fig. 9. Comparison between the observed and predicted 
amplitudes and variations in the radial- velocity, Vhigh, 
Viow, and BIS. The best fits to the observed data are 
shown (left), as well as the amplitude of the variations 
(A); the fitted curves were superposed to the simulated 
points for comparison. For this simulation, the parameters 
of the secondary CCF were 72 = — 2.1kms~^, A2 =0.14, 
(72 =8.0kms^^ and K2 =4.0kms~^. The results show the 
fits to be very good 



of the primary CCF and the 7-velocity of the secondary 
CCF (72, varied between —0.5 and — 3.0kms~^). In all 
simulations, the width and depth of the primary dip (cti 
and Al) were fixed at 4.36 kms~^ and 0.24, corresponding 
to the average parameters of the observed CCF of this star 
(typical parameters for a solar metallicity early K dwarf 
with low V sin i - see Appendix) . 

3.3. Qualitative results 

A good example of a successful simulation can be seen 
in Figs. ^ and |l^. In the former we compare the ob- 
served (left) and simulated (right) radial- velocities, Vhigh, 
Viow, and BIS for a simulation where the parameters 
of the secondary dip are 72 = — 2.1kms^^, A2 =0.14, 
(72 =8.0 km s^^ and K2 =4.0kms~^. The match is ex- 



This value is based, for simplicity, on the magnitude differ- 
ence between the two objects. 

In any case, the real shape of the bisector of HD 41004 A is 
not known. 



tremely good. In Fig. 10 we compare for the same simu- 
lation the slope of the relation BIS vs. Vr and the actual 
shape of the bisector for the maximum and minimum ve- 
locities observed and simulated. 

We note that unfortunately the comparison of the bi- 
sector shapes (Fig.|l^, right panels) have a bit of an un- 
certainty, given that we do not know the real shape of 
the bisector of the CCF of HD 41004 A. This seems to be 
particularly true for the bisector point nearest to contin- 
uum (see Fig.|^). We have thus introduced in Fig. 10 (lower 
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Fig. 10. Left: Comparison between the observed and sim- 
ulated slope of the relation between radial-velocity and 
BIS. Right: Comparison between the shape of the bisec- 
tor of the CCF for the maximum and minimum velocities. 
For the observed case, an average of the bisectors for the 
measurements corresponding to the maximum/minimum 
of the velocity was considered. The vertical lines are drawn 
just to facilitate a comparison. The error-bars drawn rep- 
resent the uncertainties in the observed bisector, as well 
as the errors due to the fact that the real "zero" bisector is 
not vertical as considered in the simulations, but slightly 
"C" shaped (see text for more details). The parameters of 
the secondary CCF are the same as in Fig.||. The velocity 
scale in these plots is similar to the one of Fig. ^ 

right panel) a tentative error bar in this point, representa- 
tive of the typical shift observed for the case of HD 13445 
- Fig. ^. If in one hand this particular point is not used 
to compute the BIS, it will change the best general visual 
bisector shapes. 

The important result at this point is that the plots 
clearly show that with the good set of parameters we can 
perfectly fit the observations. But to which extent can we 
constrain the results in order to characterize the properties 
of the CCF of HD 41004 B? In particular, it would be very 
interesting to obtain the values of a2 and K2 , that would 
permit us to derive the projected rotational velocity of 
HD 41004 B and the minimum mass for the companion 
orbiting it. 

3.4. Quantifying the parameters 

The results of the simulations described above have shown 
that there are some combinations of parameters for the 
secondary CCF that can match the observed amplitudes 
(in radial- velocity, BIS, Vhigh and Viow), and the slope 



Fig. 11. Slope of the relation between BIS and Vr against 
(72 for all the simulations run. We can clearly see the trend 
of increasing slope with increasing a2 ■ The larger dots rep- 
resent only the simulations whose results are within the 
observational constraints (see text for further details). The 
line denotes the observed slope. We clearly see that CT2 
must have values close to 8kms^^ 



of the BIS vs. Vr relation. To constrain the models, we 
have first selected all the simulations reproducing these 
quantities within 2-a (corresponding to the uncertainty in 
the observed values, i.e. 4ms~^ for the amplitudes and 
0.06 for the slope). 

We have then further compared the shape and 
"stretching" of the bisectors in order to better constraint 
the models. For this latter quantity, we have defined the 
difference in velocity between the upper bisector point 
(nearest to continuum, lets call it V{A) and V{A') in 
two bisectors corresponding to extreme velocity cases) and 
lower bisector point {V{B) or V{B')) - see Fig.|l^. The 
sum\V{A)-V{B)\ + \V{A')-V{B')l for example, is quite 
independent of the initial bisector shape, since it measures 
a differential displacement]^. Based on Fig. |^ we have con- 
sidered that we can estimate the observational values of 
\V{A) ~ V{B)\ and \V{A') ~ V{B')\ with a precision of 
50 ms^^ (to take into account the bisector "C" shape), 
and that \V{A) - V{B)\ + \V{A') - V{B')\ can be ob- 
tained with an uncertainty of 30 m s^"'^, a value that seems 
reasonably conservative. But we note again that the abso- 
lute shape of the bisectors cannot be used as a comparison, 
since the real shape of the bisector for HD 41004 A is not 
known. 



This is true since the bisector is probably not strongly "C' 
shaped for a star of this spectral type. 
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Fig. 12. Plot of the observed (solid lines) and the sim- 
ulated (dotted lines) bisectors illustrating the effect of 
varying the amplitude K2 in the simulations from 2.5 to 
lOkms^^. All the other parameters are as in Fig.|| and 

Hol. 



Although there is some degeneracy in the final results, 
some constraints can be set on the parameters of the sec- 
ondary CCF. First, we have verified that the slope of the 
plot in Fig. |l^ is most of all dependent on the width of 
the small dip. In all the simulations satisfying the cri- 
teria described above, ct2 had values between 7.4kms~^ 
and 8.6kms~^. This result, perfectly seen in Fig.|ll|, per- 
mits to determine a2 with a good precision: we estimate 
(72=8±1 kms~^. Actually, we can interpret this depen- 
dence in light of the fact that the width of the smaller 
CCF will be the most important parameter controlling 
the extent to which a change in relative position between 
the two functions will induce a change in the measured 
radial-velocity (by fitting a Gaussian function to the final 
CCF) and in the BIS=Vhigh-Viow. 

The knowledge of the width of the CCF of HD 41004 B 
means that we can estimate the projected rotational veloc- 
ity of this star. Considering that the CCF of a non- rotating 
M dwarf with solar-metallicity is 



j5kms ^ (Delfosse 
et al. |1998a| )Fl, using the calibration presented in the 
Appendix we obtain a value of ~12kms~^ for the v sini 
of HD 41004 B. As we will see below, this value will per- 
mit to set constraints on the mass of the companion to 
this star (responsible for the 1.3-day period observed in 
radial- velocity) . 



For the other parameters {K2, A2, and 72), the con- 
straints are not so strong. For example, there seems to 
be a clear relation that permits to have good models in- 
creasing A2 and decreasing if 2- It would be very inter- 
esting, however, to derive K2, since its value could give 
us an estimate for the mass of the companion orbiting 
HD 41004 B|^. First, we have verified that it cannot have 
a very high value, or else we would find a flattening in one 
(or both) extremes of the radial- velocity curve (as already 
a bit noticed in the minimum of the the simulated veloci- 
ties in Fig.^, since the smaller CCF would be getting out 
of the bigger one. Furthermore, if K2 were very high, the 
bisector shapes seen in Fig. [lo| would be quite symmetrical 
(see Fig.|l^). Based on the simulations satisfying all the 
criteria described above we can set an upper limit for K2 
around 5.2kms~^ (all the bona-fide simulations gave K2 
between 2.8knis^^ and 5.2kms^^). From now on we will 
use the "worst case" value of 5.2kms^-'^, i.e., the value 
that will correspond to the higher mass for the compan- 
ion to HD 41004 B. As we will see below, this value will 
permit to estimate a tentative upper limit for the mass of 
the companion to HD41004B. 

We could, in principle, try to see the variation in the 
observed CCF of HD 41004 AB due to the fact that it is 
constituted of the sum of two CCF's, one of which changes 
its position. To quantify the variation expected, we have 
subtracted two simulated "composed" CCF's correspond- 
ing to the extreme velocity cases of the position of the 
smaller one. For the parameters determined from our sim- 
ulation, we can show that the difference is of the order of 
0.3% (maximum). This small variation, within the obser- 
vational errors of our measurements, is simply the result 
of the fact that the CCF of HD 41004 B is very large and 
shallow. 

It is interesting to say that an analysis of the other 
parameters (width and depth) of the observed cross- 
correlation function of HD 41004 AB has revealed no spe- 
cial periodicity nor scatter. Such variations could be ex- 
pected for a typical spectroscopic binary as a consequence 
of the variation of the relative position of the CCF's of the 
two stars. In agreement with the observations, our simu- 
lations show that these parameters do not change signif- 
icantly. This is basically due to a combination of having 
a-2>>o'i and A2«Ai. 



4. A brown dwarf around HD 41004 B? 

As shown in the previous section, the best way to explain 
the observed radial-velocity signal and CCF-bisector vari- 
ations is to consider that HD 41004 B has a companion in a 
~1.3 day period orbit; considering a stellar mass of 0.4 Mq 
(typical for a M2 dwarf^ this value corresponds to a sepa- 
ration of only 0.017 AUp^. Taking the maximum amplitude 



Such a value was obtained for the ELODIE spectrograph; 
with CORALIE, the value would be a bit smaller, but given 
the large obtained (T2, the precise knowledge of this value is 
not crucial. 



In this sense 72 and A2 are not particularly important to 
constraint. 

I.e., ~7.5 stellar radius, considering that HD 41004 B has a 
radius of 0.5 Rq, typical for a M2 dwarf. 
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K2 < 5.2kms~^ deduced from the simulations above, we 
can compute an upper limit for the minimum mass of the 
companion. The results give a value of ^16 Mjup, strongly 
suggesting that the companion is in the typical brown- 
dwarf regime. We note, however, that according to our 
models, the mass could even be loweif^, and we cannot 
completely exclude that it is in the planetary regime. At 
this moment we prefer to remain cautious on this point, 
and we will thus consider this upper limit in the rest of the 
discussion; in this case, it is interesting to further discuss 
the implications of the discovery. 

For these calculations we have considered a circular 
orbit. This is not only a direct result of the observed and 
simulated radial- velocity variations, but it is in fact ex- 
pected for such a system. Actually, we can estimate the 
circularization timescale due to tidal dissipation in the 
companion's convective envelope (Rasio et al. 1996 ) - as 
done for HD 162020 by |Udry et al. (2002| ). TheTalue ob- 
tained is of the order of 10^ yr, i.e., much shorter than the 
stellar age. In other words, we expect the system to be 
circularized^. 

The presence of a third body, in this case HD 41004 A, 
may induce eccentricity growth in the system composed 



of HD 41004 B-hbrown dwarf ( |Mazeh fc Shaham 1979| ). If 
the period of this eccentricity modulation is shorter than 
the tidal circularization period, we could eventually ex- 
pect that the semi-major axis of the short period system 
would decrease, simply because of the tidal dissipation 
occurring during the circularization of the orbit. To try 
to verify if the system could have survived to this effect, 
we have computed, using Eq. 3 in Mazeh & Shaham, the 
eccentricity modulation period expected for this system. 
The result gives a value around l.STO^yr. This is both 
longer than the circularization timescale and of the order 
of the derived stellar age, and thus we can expect that the 
system constituted of HD 41004 B-|-brown dwarf is stable. 
We further note that the timescale for orbital decay due 
to tidal dissipation in the stellar convective envelope is of 
the order of lO^'^yr ( Rasio et al. 1996 ). 

We can also determine the stellar synchronization 
timescale for this system (HD 41004 B+brown dwarf). 
Using the equations in |Udry ct al. (2002D 0, we have 



obtained a value of 10^^** yr. It is thus very likely that 
this system is synchronized. We can use this fact to esti- 
mate the "true" rotational velocity of HD 41004 B, since 



^"^ Some models that fit the observations have K2 as low as 
2.8kms~^, which would imply a value for the minimum mass 
of ~8 Mjup. 

We note that the circularization timesc ale due to tidal dis- 
sipation in the stellar convective envelope ( Zahn 1989|) is much 
longer (~10^° yr)^ 



Taken from (Zahn 1989, 1992). For our case, we have con- 
and k for a 0.6 Mq star, the lower stellar 



sidered the A2, tj^ ^ 

mass tabled in ^ahn (1994 ). Since the most important factor 
determining the synchronization rate is the ratio between the 
distance between the two objects and the radius of the star, 
and that these values do not change very steeply with stellar 
mass, we believe this is a good approximation. 



the synchronization implies that the rotational period of 
the star is the same as the orbital period of the companion 
(1.3 day). Taking a stellar radius typical for a M2 dwarf 
(^0.5 Rq), we obtain a value of '^20kms^^. 

As we saw above, the width of the CCF of HD 41004 B 
can be fairly well constrained by our simulations, per- 
mitting to make a good estimate for the projected ro- 
tational velocity v sinz~12.0kms~^. This value implies a 
sin?~0.6. Within these conditions, we can deduce that 
the companion to HD 41004 B has a mass lower than 
^25Mjup, placing it in the usually considered brown- 
dwarf regime. 

At a separation of only ^3.7 Solar radius, we can cal- 
culate that the Roche lobe of th e companion is not filled. 
Using Eq. 2 in Egglcton (1983 ), we obtained a value of 
RRL=0.18a, where a is the semi-major axis of the orbit. 
With a=3.7RQ, we obtain Rrl'~^0.7Rq. This value is 
a factor of 7 larger than the typical radius of a 25Mjup 
brown dwarfQ. No mass transfer is thus expected to occur. 



5. Concluding remarks 

We have presented the case of HD 41004 AB, a system 
composed of a KOV star and a 3.7 magnitudes fainter M2 
dwarf, separated by 0.5". Radial-velocity measurements 
derived from CORALIE blended spectra of the two stars 
have unveiled a radial-velocity variation with a period of 
~1.3days and a small amplitude ('^SOms"^), compatible 
with the expected signal due to the presence of a planetary 
companion to HD 41004 A. However, as we have seen, the 
combined radial-velocity, photometry and bisector analy- 
sis suggest that the best explanation for the observations 
is that the fainter HD 41004 B has a brown-dwarf compan- 
ion. In this scenario, the observed low amplitude radial- 
velocity variation is due to the measure of a variation in 
the line profiles which is induced by the change of the 
relative position of the spectra corresponding to the two 
stellar components HD 41004 A and B. 

If confirmed, the present discovery represents the 
first detection of a short-period brown-dwarf companion 
around a M2 dwarf. In particular, its estimated upper 
limit mass (~25Mjup) puts it in the middle of the so- 
called brown-dwarf desert, a mass region (between ~20 
and 40Mjup) for which (alm ost) no short period co mpan- 



ions to solar-type F, G, K ( Halbwachs ct al. 2000 Udry 



3t al 200l| ; [Jorissen et al. 20011 ) and M ( pVIarcy fc Benitz 



1989) dwarfs were found. The fact that HD 41004 B is a M 



dwarf makes us speculate that the formation of such sys- 
tems is more likely for lower mass primaries, i.e. systems 



having mass ratios closer to unity (Duquennoy & Mayor 



1991) 



Armitage & hionneli {'ZWZ) proposed a model to ex- 
plain the existence of the brown-dwarf desert. One inter- 
esting feature of their model is that it predicts that no 
brown-dwarf desert should be observed for the compan- 



At an age of ~lGyr, a brown dwa rf with a mass aroun d 
0.03 Mq has a radius around 0.1 Rq - Chabrier et al. (200C). 
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ions around the lowest mass dwarfs; they set an upper 
limit of 0. 1-0.2 M0. Besides the fact that HD 41004 B is 
slightly more massive than this limit, the present case is 
interestingly similar to the predictions. 

Recently, Zucker fc Mazeh (2002 ) discussed an inter- 
esting correlation between the planetary mass and the or- 
bital period. Their analysis strongly suggests that for sin- 
gle stellar systems, there is a lack of "high" mass planetary 
companions in short period orbits. On the other hand, 
Zuckcr & Mazeh have found that for stars in multiple sys- 
tems, this correlation is no longer present. In particular, 
they have pointed out that in these latter cases, there 
seems to be a negative correlation between "planetary" 
mass and orbital period (see their Fig. 3). It is very inter- 
esting to see that the companion to HD 41004 B perfectly 
fits this trend. 

A few works have tried to study the formation of plan- 
ets (or low ma ss objects) around stars in multiple systems. 
Nelson (2000 ) showed that the formation of a planet in a 
disk is unlikely for equal mass binary systems with a sepa - 
ration lower than ^50 AU. On the other hand. Boss (1998 ) 
suggested that the influence of a companion might trig- 
ger "planetary" formation by disk instability. Although no 
strong conclusions are possible at this moment, the fact 
that HD 41004 B is in a double system with a separation 
that can be as low as ^20 AU is very interesting from the 
point of view of the formation of its companion. 

Given the uncertainties in the shape of the observed bi- 
sector and in the models, we do not pretend to have a pre- 
cise mass determination for the companion to HD 41004 B. 
Although we believe we have obtained a good estimate, 
one of the main goals of this paper was to illustrate the 
importance of combining the radial- velocity data with the 
bisector analysis. In other words, one important lesson 
to be taken from the presented results is that the bisector 
analysis was crucial to correctly interpret the observations. 

It is important to caution that this kind of situation 
can also happen for long period systems. For those, the 
signal is very unlikely to have an intrinsic stellar activity 
origin, but we cannot exclude that it might originate from 
the presence of a "wobbling" stellar companion. It is very 
easy to find a situation where an undetected companion, a 
few tens of an arcsec distant, can be "contaminating" our 
analysis. Together with a good knowledge of the target 
star environment, in such a situation the bisector analysis 
seems to provide an unique tool to point out the exact 
origin of the radial-velocity variations. 

In this sense we have analyzed (or re-analyzed) the 
bisectors for all the stars with planets discovered in the 
context of the CORALIE planet search programme^. We 
did not find any significant correlation between Vr and 
BIS. We can thus remain confident that the presence of 
a planetary companion is the best way of explaining the 
radial- velocity variations in these systems. 

Besides the bisector analysis, one other way of test- 
ing a radial-velocity variation for such cases may in- 
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Fig. 13. Phase folded radial-velocity measurements of 
HD 41004 AB obtained using two diflFerent CCF masks. 
As expected, the amplitudes K obtained vary from mask 
to mask (see text). The vertical scales of the two plots is 
the same to facilitate a comparison. 

volve the measurement of the radial-velocity using dif- 
ferent sets of lines, or different spectral regions. If the 
radial-velocity variation is due to the presence of a planet, 
we can expect that every spectral region/line will give 
us about the same velocity amplitude (but not necessar- 
ily the same 7- velocity). The fact that for HD 41004 AB 
the amplitude in radial velocity was different using the 
cross-correlat ion mask constructe d specially for the bisec- 
tor analysis ( Queloz et al. 2001b| ) from the one obt aine d 
using the "classical" mask is very telling (see Sect. 2.3). 
In Fig. |l^ we can further see two phase folded diagrams 
for the radial- velocities of HD 41004 AB obtained using 
two different CCF masks. As we can see, if we use a mask 
specially constructed for the radial- velocity determination 
of M4 dwarfs ( [Delfosse et al. 1998b ), i.e. a spectral type 
close to the one of HD 41004 B, we obtain a much higher 
amplitude in radial-velocity than for the case of using a 
mask constructed for KO dwarfs. This difference is ex- 
pected since in such a case, the difference between the 
CCF's of HD 41004 A and HD 41004 B is much smaller (the 
M dwarf spectrum is enhanced relatively to the K dwarf), 
being thus the influence of the "small" CCF stronger. A 
similar situation (although not as strong) is seen when 
comparing the amplitudes obtained using the KO and the 
FO masks. 

This kind of analysis can, in principle, serve as a test, 
if the study of the bisector is not possible. But we note 
that for radial- velocity variations induced by the presence 
of dark spots in the stellar photosphere, i.e. where the 
same phenomenon is affecting all spectral lines in about 
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the same way, we do not expect a strong difference be- 
tween the radial-velocity amplitudes obtained using dif- 
ferent spectral lines. 

Finally, it is interesting to say a few words about other 
possible ways of confirming the current detection. Given 
the short period of the brown dwarf around HD41004B, 
the probability that we are able to observe a transit is 
quite high. In a first approximation, the magnitude vari- 
ation expected in such a transit is around 10^^ (for the 
A-fB system), a value that does not seem too low. But a 
simpler way of confirming this case would pass by doing 
high-resolution spectroscopy (and velocity measurements) 
of HD 41004 B. However, this is not an easy task, since the 
two components of HD 41004 AB are separated by only 
0.5". Unfortunately, there are no available high-resolution 
spectrographs attached to an Adaptative Optics system 
in the southern hemisphere capable of accomplishing this 
task. The use of the Hubble Space Telescope (HST) might 
represent a solution. Else, the solution may pass by using 
high-resolution ncar-IR spectroscopy, since at those wave- 
lengths the flux of HD 41004 B is much closer to the one 
from HD 41004 A. 



0.95 - 



0.9 - 







n 


hi KOV 




J:j [Fe/H] = -1.8 


■ 1 


Sun (G2V) 




[Fe/H]«0.0 




i KOV 


I.I. 


J [Fe/H] = 0.0 - 
1 . 1 



-40 



(V - Vr) [km/s] 



-20 20 
(V - Vr) [km/s] 



40 



Fig. A.l. Left panels: plot of the CORALIE CCF for two 
stars with different projected rotational velocities. This 
figure illustrates the effect of the v sini on the CCF; right 
panel: plot illustrating the variation of the surface of the 
CCF due to changes in the metallicity and temperature 



can be used to derive two important quantities from the 
CORALIE CCF: the v sini and the stellar metallicityQ. 



Appendix A: Determination of v sini and [Fe/H] 
from the Cross-Correlation Function 

One interesting property of the Cross-Correlation 
Function (CCF) is that any stellar "phenomenon" capa- 
ble of infiucncing the lines included in the correlation mask 
(e.g. the line profiles) will be reflected on the CCF itself 
(i.e. on its shape, depth and width). In other words, the 
CCF represents an "average" spectral line amidst those 
used in the correlation mask (Mayor 1985). The profile 



and intensity of the CCF is thus a convolution between the 
profile due to the intrinsic stellar atmospheric parameters 
affecting the lines in the mask (like the global abundance, 
thermal broadening, pressure broadening, or microturbu- 
lence), and the macroscopic macroturbulence and rota- 
tional velocity of the star ( Gray 1992 ), without forgetting 
the instrumental profile, characteristic of a given instru- 
ment. In the case of CORALIE, given that the correlation 
mask is constructed using mostly weak neutral lines in 
the spectrum of a standard star (a KO dwarf), the prop- 
erties of the resulting CCF are basically dependent on the 
physical quantities controlling the properties of a "typi- 
cal" neutral weak metal line. For a detailed description of 



the Cross-Correlation technique we refer to Baranne et al 
p996| ) and |Pepc et al. (200^ ). 



From a practical point of view, this give us the pos- 
sibility to obtain stellar quantities that are reflected on 
the spectral lines directly by analyzing the parameters of 
the CCF. We can even hope to build specific masks for 
different kinds of lines, each sensitive to a different stel- 
lar parameter. Let us see in the next sections how this 



A.l. A calibration of the projected rotational velocity 

V sini 

The Gaussian width of a weak spectral line of a "non- 
rotator" depends basically on the spectral type and lu- 
minosity class. This is mostly related to the fact that 
temperature and surface gravity are the main variables 
controlling both the line strength and broadening, due to 
the temperature and pressure effects^. These variables 
actually indirectly control other broadening parameters 
like the macro- and micro-turbulence ( Gray 1992 ). At the 
end, the only broadening parameter left that is definitely 
independent of the intrinsic stellar atmosphere properties 
is the projected rotational velocity {v sini). 

Actually, considering that the instrumental profile is 
constant for each star and does not change in time for 
a given spectral type and luminosity class (which is a 
reasonable thing to do for such stable instruments like 
CORALIE) we can expect that every star will have ap- 
proximately the same "zero" CCF width - hereafter called 
(To ^ and for a constant temperature and surface gravity, 
the "excess" width that the lines (i.e. the CCF) might 
present is then mostly due to the effect of the rotational 
velocity. 



For more details, and a broader discussion about the errors 



for th e two calibrations presented below, we refer to Santos 



(2002). 

Two other variables, the abundance and the magnetic 
field, will be discussed below. Note also that we will "forget" 
about natural broadening, since it is constant for each line. 
Furthermore, we consider that the limb-darkening properties 
are the same for a given temperature and surface gravity. 
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Fig. A.2. Plot of the Gaussian width of the CORALIE 
CCF (expressed in kms~^) vs. B — V. The lower enve- 
lope (solid line) determines the locus of the stars with 
zero rotational velocity. The two dashed lines represent 
the uncertainty of 0.025 km s^^ in the position of the lower 
envelope of the points 



There are, however, two other physical quantities that 
can (in principle) change significantly the value of ag: the 
metallicity and the magnetic field. In fact, the higher the 
metal content of a star, the higher the number of satu- 
rated lines: this saturation will cause a slight increa,se in 



(To- Tn lour case this effect is not very strong, since the 
cross-correlation masks are based mainly on weak lines. 
As for the magnetic field, it is known that changing its 
intensity will change the Zceman spliting of the spectral 
lines, t hus broadening their profiles. This fact can even be 



that the CCF of a rotating star can be approximated by 
a Gaussian. As shown by Queloz et al. (199S) this is true 
up to velocities of ~20kms^^, the regime we are mostly 
interested in. 

Using the same idea, we have used the CORALIE CCF 
to calibrate a relation between the a obtained with this 
instrument and the v sini for solar type dwarfs. There are 
basically two things to do. The first is the determination 
of (Tq. This quantity can be obtained simply by adjusting 
the lower envelope of the distribution of points in a plot of 
(T vs. B — V (this latter quantity is quite well related with 
temperature)^ In this envelope we will, in principle, find 
the non-rotator stars. We note that most of the stars ob- 
served with CORALIE are dwarfs, and there are probably 
no objects with luminosity classes lower than IV. Such a 



plot is seen in Fig. A. 2, and the adjusted lower envelope is 
well described by: 



(A.2) 



o-o = 6.603- 6.357(5- V^) 

+5.533 {B - Vf - 1.454 {B - Vf 



It is interesting to say a few words about Eqn. [A.2| . As 
can be seen from Fig. A.2, the fitted function decreases up 
to values of i? — V^^O.9 (where it has a minimum) and then 
increases again towards higher B — V values. In the low 
B — V regime, this decrease might simply be explained by 
a decrease in the macroturbulent dispersion ( pray 1992 ). 
On the other side, for high B — V values, the increase in ctq 
is probably due to the increase in the damping constant 
due to van der Waals broadening, that can be shown to 
increase with decreasing temperature. Furthermore, up to 
a spectral type of K5, neutral metallic lines increase in 
strength, thus increasing the number of saturated lines 



(for a constant abundance). As discussed in Benz & Mayor 



(1984 ), the increase of the magnetic field as a funcion of 



used td obtain an estimation of the stellar magnetic field 
(e.g. Borra et al. 1984 ; Queloz et al. 199(: ). However, it 
has been shown by Benz fc Mayor (1984 ) that the line 



increasing spectral type might also be responsible for the 
increase in ctq after a given spectral type. Together with 
these physical processes, variations resulting from the way 
the CCF mask is optimized (Baranne et al. 197' 



Benz 



Mayor 1984) might also introduce some trendm. It is. 



broadening due to the effect of the magnetic field is only 
important for late K and M stars. Since this effect would 
be quite complicated to take into account, and it is prob- 
ably not important for most of the stars in the CORALIE 
survey, we will not consider it in the rest of this discussion. 

Benz & Mayor (1984) have successfully shown that the 
width of the CORAVEL CCF is well correlated with the 
V sini of a star. These two variables can be related by: 



however, difficult to quantify exactly what is the net effect 
of all these processes. 

It is important to mention that after all the metallic- 
ity does not seem to play an important role in the de- 
termination of g p . A few tests were done, separating the 
stars in Fig. |A.2| into different metallicity bins. The result 
showed that the effect is very small, and difficult to quan- 
tify amidst the noise in the lower envelope of the points. 
The interpretation of this is, as mentioned above, proba- 
bly quite simple: since we are dealing with weak lines in 



A 



(A.l) 



where cr represents the measured Gaussian width of the 
CCF, ctq the value of the expected a for a "non-rotator" of 
a given spectral type and luminosity, and A is a constant 
relating the "cr-excess" to the actual projected rotational 
velocity of the star (u sini). We are, of course, considering 



Since the B — is a widespread index available for all 
the stars in the CORALIE sample, it is a particularly suitable 
variable to use. 

We expect a to be slightly different in the red and in the 
blue part of the spectrum since the "lines" in the mask are 
in average slightly larger in the red; consequently, this might 
introduce a small increase of the CCF with increasing spectral 
type. 
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Fig. A.3. Plot of the CORALIE CCF surface Wf^t as a 
function oi B — V. Wfu, expressed in kms"^, is defined 
as •\/2 TT • a fit ■ ^fit, where afu and A fa are the Gaussian 
width (in km s^^) and depth of the measured CCF, respec- 
tively. The lines represent regions with [Fe/H]=constant. 
The size of the points is proportional to the spectroscopi- 
cally determined metallicity. The spectroscopic metallici- 
ties were taken from Santos et al. (2001a . b[) 



the linear part of the curve of growth, the saturation effect 
is not expected to be very strong. 



The other variable in equation AT we need to deter- 
mine is A, i.e., the constant relating the v sinz to the 
"excess" width of the CCF. This can be done in the same 



way as in Queloz et al. (1998 ): by measuring the variation 
of a as we convolve the CCF of non -rotating stars with 
a given rotational profile ( Gray 1992 ). For CORALIE we 
find A=1.9±0.1, close to the value obtained for ELODIE 



( Queloz et al. 1998 ). In any case, the final value of v sinz 
is not strongly dependent on errors in A. 

The technique described above represents a very sim- 
ple way of obtaining precise projected rotational velocities 
for dwarfs, simply as a by-product of the precise radial- 
velocity measurements. In particular, the calibration per- 
mits us to obtain quite easily values for the v sini for all 
the stars in the CORALIE planet search programme. 



A. 2. A calibration of the metallicity 

Another variable that can easily be obtained from the 
CCF of a star is its metallicity. This has been firstly done 



by [Mayor (1980[ ) for the C ORAVEL C CF, a cahbration 
that was later improved by Pont (1997). 

As mentioned above, the mask used for CORALIE is 
mainly built out of a set of weak neutral spectral lines 



Fig. A. 4. Upper panel: Comparison between the spectro- 
scopic and calibrated [Fe/H]. Lower panel: Residuals of 
the calibration as a function oi B — V; no dependence 
is found. The spectroscopic metallicities were taken from 
[Santos et al. (200l4 [b|) 

(from a template KO spectra). Since most of them are 
iron lines (the main contributor for the line opacity in 
solar type stars), we can expect that the surface (i.e. the 
Equivalent Width) of the CCF is well related to the [Fe/H] 
of a star. 

Furthermore, the fact that we are dealing with weak 
neutral lines in solar-type stars implies that the surface 
of the CCF will be basically independent on the micro- 
turbulence and the surface gravity. These would have an 
important effect if we were not dealing with lines that 
are in the linear part of the curve of growth, since the 
broadening produced by these effects would act to de- 
saturate the lines, increasing thus their equivalent widths. 
Furthermore, it can be shown that weak lines of an ele- 
ment for which most of it is in the next ionization state 
(most of the iron in a solar- type star is in the form of Fe ll) 
are quite insensitive to pressure changes. The surface of 
the CCF will thus mainly depend on the temperature (well 
correlated with B — V) and abundance (see Fig A.l). 

To calibrate a relation between the surface of the 
CORALIE CCF (hereafter Wm) and the iron abundance 
expressed by [Fe/H], we have used a set of stars for which 
we had obtained precise [Fe/H] values using a detailed 
spectroscopic analysis ( pantos et al. 2001a ^). The fit, 
shown in Fig. A. 3 gives: 



[Fe/H] 



2.573 + 4.587 \ogWcor 

-8.142 {B-V) + 3.583 {B - V)"^ 



(A.3) 



a calibration valid for 0.52 < B-V < 1.1, 1.18 < Wm < 
3.59 and -0.55 < [Fe/H] < 0.42. As expected, the figure 
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shows that for a given metalhcity, the equivalent width of 
the spectral lines (i.e. the surface of the CCF) increases 
with decreasing temperature. This is the "normal" be- 
haviour of a weak metallic line in this temperature regime 
(up to spectral types around K5, B — V^l.2). 

A comparison between the metallicities obtained from 
the CCF surface and a spectroscopic analysis shows that 
the fit has a remarkable small dispersion (0.05 dex - 
Fig. A.4), similar to the uncertainties in the spectroscopic 
determinations. Furthermore, there seems to be no spe- 
cial systematics with temperature. We can thus trust on 
this calibration to obtain precise values of [Fe/H] for the 
stars in the CORALIE planet search programme without 
any long and arduous spectroscopic analysis. We note that 
the metallicities obtained this way are spectroscopic de- 
terminations, since we are using spectral line information. 
Although we cannot exclude that for a particular star the 
result might not be as precise as a detailed spectroscopic 
analysis determination, at least in statistical studies this 
is certainly an extremely accurate technique. 
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